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We review the current status of the study of parity and time invariance violation in 
atoms, nuclei and molecules. We focus on parity non-conservation in cesium and three 
of the most promising areas of research: (i) parity non-conservation in a chain of isotopes, 
(ii) search for nuclear anapole moments, and (iii) search for permanent electric dipole 
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EDM or nuclear T, P-odd moments such as a nuclear EDM or nuclear Schiff moment. 
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1. Introduction 

The study of parity and time invariance violation at low-energy in atoms, molecules 
and nuclei is a relatively inexpensive alternative to high-energy searches for new 
physics beyond the standard model (see, e.g. reviews [TJ [2]) . Accurate measurements 
of the parity non-conservation (PNC) in atoms is one of the most promising ways of 
exploring this path. It culminated in a very precise measurement of the PNC in the 
cesium atomP There was even an indication that this measurement showed some 
disagreement with the standard model and might indeed lead to new physicsP It 
took considerable efforts of several groups of theorists to improve the interpretation 
of the measurement and resolve the disagreement in favor of the standard model. 
The disagreement was removed when the BreitP and quantum electrodynamic cor- 
rections^ were included and the accuracy of the treatment of atomic correlations 
was improvedP^ 

The latter worlP claimed a very small theoretical uncertainty of just 0.26%. 
Combined with the PNC measuremenlP it leads to perfect agreement between the 
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measured weak nuclear charge of the 133 Cs nucleus and the value predicted by the 
standard model. It also puts strong constraints on new physics beyond it. 

The PNC amplitude calculated in Ref. |8] is about 0.9% smaller than that of 
earlier calculations EEE] The difference is significant given that all earlier works 
agree very well with each other and theoretical uncertainty in them is estimated to 
be 0.4 - 0.5%P^ The effect of this difference on the agreement with the standard 
model is not large, however, there are important implications for the constraints 
on new physics beyond it. The difference was attributed in Ref. [8] to the role of 
high-order correlations. However, recently in Ref. [12] it was found that the authors 
of Ref. 8 overlooked some important correlation contributions to the PNC ampli- 
tude. The uncertainty of 0.26% claimed in Ref. |8] is based mostly on analysis of 
the dominating contributions. The uncertainty due to the correlation corrections 
in other terms has been underestimated. The result of Ref. |8| when corrected as 
suggested in Ref. [TH is in perfect agreement with the previous calculations EEHU] 
PNC in Cs is discussed in more detail in the next section. 

It is unlikely that any new measurements of the PNC in atoms can compete 
with the cesium experiment in accuracy of the interpretation. There are conflict- 
ing tendencies in the accuracy of the measurements and the calculations. On the 
one hand, it is easier to get highly accurate measurements for heavy atomic sys- 
tems where the PNC effect is larger. On the other hand, advanced computational 
techniques favor light alkaline atoms. The best accuracy for both calculations and 
measurements have been achieved for cesium atom: 0.35% for the measurements^ 
and 0.4 - 0.5% for the calculations EEK121 A higher accuracy of the interpretation is 
possible for the rubidium atom, which is a lighter analog of cesium. The rubidium 
atom has been considered so far for a different type of PNC effect: the nuclear spin- 
dependent PNC in the ground-state hyperfme transition, which is mostly due to 
the nuclear anapole momenlP^ (see section 4) . For the spin-dependent PNC accu- 
rate calculations are also possible, but less important than for measurements of the 
weak nuclear charge. It is unclear whether the weak nuclear charge of rubidium can 
be measured to high accuracy. A heavier analog of cesium, the francium atom, is 
considered for both types of PNC measurements, the effect of weak nuclear charge 
and the anapole moment ESEH The accuracy of calculations for francium could be, 
at best, as good as for cesium, but may be slightly lower due to the larger relativis- 
tic effects which increase the Breit, quantum electrodynamics (QED) and neutron 
distribution ("neutron skin") corrections. There are suggestions to measure PNC in 
s-d transitions of Ba+G^ 1 and Ra+GsE 3 

ions which have an electronic structure 
similar to that of the cesium atom. Here again the accuracy of the calculations is 
unlikely to outperform cesium due to larger correlations in d-states (only s and p 
states are involved in cesium) . Note that some improvement here may be achieved 
using experimental values of the electromagnetic El p — d amplitudes in the sum 
over intermediate p states in the calculation of the s — d PNC amplitude^ . 

However, there are areas of symmetry breaking in atomic systems where we do 
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Table 1. Correlated PNC amplitude in Cs 
[in 10 — Qw/N) a.u.] in different cal- 
culations. Breit, QED and neutron skin 
corrections arc not included. 



Value Source and method 



0.908(9) CP+PTSCI, Ref.[9] 

0.909(9) CC SD, Ref.flO] 

0.905(9) MBPT with fitting, Ref.lTTI 

0.9078(45) CP+PTSCI, Ref.0 

0.8998(24) CC SDvT, Ref.® 



not need very high accuracy of atomic calculations: (i) the PNC measurements for 
a chain of isotopes; (ii) the measurements of nuclear anapole moments; and (iii) 
the measurements of the P,T-odd permanent electric dipole moments of atoms and 
molecules. Accurate numerical atomic calculations are not needed for the interpre- 
tation of the PNC measurements for a chain of isotopes since the change of the 
atomic wave function is very small and may be estimated analytically in the ratio 
of PNC effects for different isotopes. The interpretations of the anapole and electric 
dipole moment measurements do rely on atomic or molecular calculations. How- 
ever, high accuracy here is not as important as for the measurements of the weak 
nuclear charge (where we search for a small deviation from the Standard model 
predictions). 

Below we briefly review parity non-conservation in cesium, and each of the three 
subjects listed above. 

2. Parity non-conservation in cesium 

Parity non-conservation in cesium is currently the best low-energy test of the elec- 
troweak theory. It is due to the high accuracy of the measurements and its inter- 
pretation. The experimental value for the PNC amplitude is^ 



where j3 is the vector transition amplitude. The amplitude -Epnc is proportional to 
the weak nuclear charge Qw 



where K is the electron structure factor which comes from atomic calculations, N is 
the number of neutrons (N = 78 for 133 Cs). The main challenge in the calculation 
of the PNC amplitude is an accurate treatment of the inter-electron correlations. 
There are also minor contributions arising due to the effects of the Breit interac- 
tion, quantum electrodynamics and neutron skin corrections. We first consider the 
dominant term (the correlated PNC amplitude) and add minor corrections later. 

Table [1] shows the most accurate values of the correlated amplitude found in 
different calculations. The abbreviation CP+PTSCI stands for the correlation po- 
tential combined with the perturbation theory in screened Coulomb interaction, 



Im(Ep NC )/P = 1.5935 (56) mV/cm. 



(1) 



£pnc = K(-Q w /N), 



(2) 
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Table 2. All significant contributions to the Bpnc [i n 
10" 11 i(-Q w /N) a - u -l for Cs. 



Contribution 


Value 


Source 




Ref . [8] Rcf . EU 




Core (n < 6) 


-0.0020 (2) 0.0018 (8) 


Ref. H[12] 


Main (n = 6 - 9) 


0.8823 (17) 


Ref.U 


Tail (n > 9) 


0.0195 (20) 0.0238 (35) 


Ref.m[l2] 


Subtotal 


0.8998 (25) 0.9079 (40) 


Ref.l8llT2l 


Breit 


-0.0055 (1) 


Ref.|5] 


QED 


-0.0029(3) 


Ref.U 


Neutron skin 


-0.0018(5) 


Ref. [5] 


Total 


0.8906(26) 0.8977 (40) 


Ref. El [12] 



CC SD stands for the coupled cluster with single and double excitations, CC SDvT 
means the coupled cluster with single, double and valence triple excitations. One 
can see that all results except those of Ref.[5]agree with each other within 0.4%. The 
result of Ref. [5] is 0.9% smaller than another most accurate result of Ref. [7J Tabled] 
shows the reason for this difference. The authors of Ref. [5] used the sum-over-states 
approach in which the PNC amplitude is expressed as a sum over complete set of 
single-electron states, including states in the core and states above the core. Con- 
tribution of each state is calculated independently of others and the use of the same 
approximation for all terms is practically impossible. A very sophisticated coupled 
cluster method was used to calculate the main term (contributions of the four low- 
est states above the core). Significantly less accurate approximation was used for 
the tail (contribution of the states above 9p) and a very simple Hartree-Fock ap- 
proximation was used for the core. The uncertainty of just 10% was assigned for 
the sum of the core and tail contributions. However, as it was shown in Ref. [12] the 
inclusion of the core polarization effect for the core contribution changes even the 
sign of the core contribution brining the total change to about 200%, far beyond 
the declared 10% uncertainty. The tail contribution is also significantly larger when 
the Brueckner-type correlations are included. When the core and the tail contribu- 
tions are corrected as suggested in Ref. [T2l the result of Ref. [8] comes into excellent 
agreement with the previous calculations (see "Subtotal" line in Table [2] and Table 
Q] the corresponding numbers are 0.9079 and 0.9078). 

Table [2] also lists all other significant contributions to the PNC amplitude in 
cesium. The resulting PNC amplitude is 

E PNC = 0.8977 (45) x IQ- Xl i{-Q w /N) . (3) 

Note that the uncertainty coming from the analysis of Ref. [12] is smaller (±0.004, 
see Table [2]) due to the small uncertainty of the main term claimed in Ref. [8] For 
a conservative estimate we use in (j3j a slightly larger uncertainty of Ref. [7] The 
central points of Ref. [7] and Ref. [12] are practically identical. 

To find the value of the weak nuclear charge of 133 Cs from (JTJ) one also needs the 
value for the vector transition amplitude /3. The most accurate value for ft comes 
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from the analysis^ of the Bennett and Wieman measurement ipSl 

/3 = 26.957 (51) a%, (4) 

where as is the Bohr radius. Comparing <j3j> , (p} and (|4]) leads to 

QwC 33 Gs) = -72.58 (29) oxpt (36) tho o ry . (5) 

This value is in a reasonable agreement with the prediction of the standard model, 
QSM = _73.23 (2j2H (see also Refs. E2 [23l Eg . If we add theoretical and experi- 
mental errors in ([5]) in quadrature, the Cs PNC result deviates from the standard 
model value by 1.4 c: 

Qw - Qm 1 = $Qw = 0.65 (46). (6) 

For small deviations from the Standard Model values we may relate this to the 
deviation in sin 2 6\y using the simple relationship SQw w — 4Z £(sin 2 6\y) which 
gives (S(sin 2 6 W ) = -0.0030 (21) and 

sin 2 6 W = 0.2356 (21) . (7) 

This is 1.4cr off the standard model value 0.2386 (1)P at near zero momentum 
transfer. 

The new physics originated through vacuum polarization of the gauge boson 
propagators, and is described by the weak isospin conserving S and isospin breaking 
T parameters^! 

Qw - Qm 1 = -0.800 S - 0.007 T. (8) 

At the lc level © leads to S = -0.81 (58). 

Finally, a positive AQw could also be indicative of an extra Z boson in the 
weak interactiorP^l 

Qw ~ Qlf w 0.4(27V + Z){M w /M Zx f. (9) 
Using © leads to M Zx > 700 GeV/c 2 . 

3. Chain of isotopes 

It is convenient in this section to present the values measured in atomic PNC- 
experiments in a form similar to (f2J) 

Epnc — kpNcQw, (10) 

where fcpjvc is an electron structure factor which comes from atomic calculations, 
and Qw is the weak nuclear charge. Very sophisticated calculations are needed for 
accurate interpretation of the measurements as has been discussed in the introduc- 
tion and previous section. An alternative approach was suggested in Ref. [25] If the 
same PNC effect is measured for at least two different isotopes of the same atom 
than the ratio 

i^PNC HW 
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of the PNC signals for the two isotopes does not depend on the electron structure 
factor. It was pointed out however in Ref. [26] that possible constraints on the new 
physics coming from the isotope ratio measurements are sensitive to the uncertain- 
ties in the neutron distribution which are sufficiently large to be a strong limitation 
factor on the value of such measurements. The problem was addressed in Ref. [27] 
and more recently in Ref. [28] The authors of Ref. [27] argue that experimental data 
on neutron distribution, such as, e.g. the data from the experiments with antipro- 
tonic atomspSl can be used to reduce the uncertainty. In the more general approach 
of Ref. 28 nuclear calculations are used to demonstrate that the neutron distribu- 
tions are correlated for different isotopes. This leads to significant cancelations of 
the relevant uncertainties in the ratio (fTT|) . 

The parameter T of the sensitivity of the ratio (fTTj) to new physics can be 
presented in the form 



where h p is the new physics coupling to protons (AQ new — Zh p + Nh n ), ho comes 
from the SM, IZo is the ratio (jTTj) assuming that each isotope has the same proton 
and neutron distribution (no neutron skin), N and N' are the numbers of neutrons 
in two isotopes, Z is the number of protons and A7V = N' — N. The constraints 
on new physics parameter h p are affected by the experimental error SlZ exp and 
uncertainties in IZo due to insufficient knowledge of neutron distributions. However, 
as is argued in Ref. [25] these errors are correlated and strongly cancel each other. 
Indeed, the distribution of the core neutrons is nearly the same in different isotopes. 
Significant changes happen in a small number of valence nucleons only. Estimations 
of Ref. 28 show that corresponding contribution to ST is in the range 10~ 3 -j- 10~ 2 
which is about an order of magnitude smaller than the uncorrelated one. In the 
end, the isotope-chain measurements may be more sensitive to new physics than 
current parity-violating electron scattering measurements^! (by a factor of 10 for 
such atoms as Cs, Ba and Dy). 

Experiments on isotope chains are in progress at Berkeley for Dy and Yb 
atomsJSS at TRIUMF for Fr atoms pi at Los Alamos for Yb+ ions, 33 and at 
Groningen (KVI) for Ra + ionsP^ There is an interesting recent suggestion to mea- 
sure the PNC in metastable Xe and Hg.^ Most of these experiments consider also 
the measurements of the nuclear anapole moment. 

4. Anapole moment 

The notion of the anapole moment was introduced by B. Ya. ZeldovichPsl Nuclear 
anapole moment (AM) is the magnetic P and C-odd, T-even nuclear moment caused 
by the P-odd weak nuclear forces. Interaction of electrons with AM magnetic field 
(which may be called the PNC hyperfine interaction) dominates the nuclear-spin- 
dependent contribution to the atomic or molecular PNC effect. 




(12) 
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First calculations of the nuclear AM and proposals for experimental measure- 
ments were presented in Refs. [3611401 Corrections to the AM interaction with elec- 
trons due to finite nuclear size were considered in Ref. [41] The status of the nuclear 
AM many-body calculations is presented in reviews [IJ [2j The authors of Refs. [38l 
[39l (see also Ref. |42| note in particular that the nuclear-spin-dependent PNC effect 
is strongly enhanced in diatomic molecules due to mixing of the close rotational 
states of opposite parity including mixing of A or f2 doublets. The PNC effects 
produced by the weak charge are not enhanced. Therefore the AM effect domi- 
nates PNC in molecules. This greatly simplifies the detection of AM in diatomic 
molecules compared to atoms. In atoms the AM effect is 50 times smaller than the 
weak charge effect; AM effect is separated as a small difference of the PNC effects in 
different hyperfine transitions. A review of the parity and time invariance violation 
in diatomic molecules (including the AM effect) can be found in Ref. [43] 

The idea of the AM contribution enhancement may be explained as follows. 
After the averaging over electron wave function the effective operator acting on 
the angular variables may contain three vectors: the direction of molecular axis N, 
the electron angular momentum J and nuclear spin I. Scalar products NI and NJ 
are T-odd and P-odd. Therefore, they are produced by the T, P-odd interactions 
discussed in the next section. P-odd, T-even operator Vp must be proportional to 
N[J x I]. It contains nuclear spin I, therefore, the weak charge does not contribute. 
The nuclear AM is directed along the nuclear spin I, therefore, it contributes to 
Vp . The matrix elements of N between molecular rotational states are well-known, 
they produce rotational electric dipole transitions in polar molecules. Therefore, Vp 
(induced by the magnetic interaction of the nuclear AM with molecular electrons) 
mixes close rotational-hyperfinc states of opposite parity. The interval between these 
levels is five orders of magnitude smaller than the interval between the opposite 
parity levels in atoms (by the factor m e /M where m e and M are the electron 
and reduced molecular masses), therefore PNC effects are five orders of magnitude 
larger. Further enhancement may be achieved by a reduction of the intervals by an 
external magnetic fieldP^ 

The effect is further enhanced for heavy molecules. It grows with the nuclear 
charge as Z 2 A 2 / 3 R(Za), where R[Za) is the relativistic factor which grows from 
R = 1 at low Z to Jf, ~ 10 for Z > 80 and the factor A 2 ! 3 comes from the 
nuclear anapole, A is the nucleon number. Good candidates for the measurements 
include the molecules and molecular ions with the projection of the electron angular 
momentum on the molecular axis fl = A z + S 2 = 1/2, i.e. in £1/2 (A 2 = 0) or H1/2 
(A z = 1) electronic ground statespSEal for example, YbF, BaF, HgF, PbF, LaO, 
LuO, LaS, LuS, BiO, BiS, YbO+, PbO+, BaO+, HgO+, etc. Molecular experiments 
are currently in progress at Yale^ and Groningen KVlpSl An interesting idea of 
studying AM contribution to the NMR spectra of chiral molecules was discussed in 
Ref. US 

Interpretation of the AM measurements requires electron structure calculations. 
A number of scmicmpirical and ah initio calculations have been performed for di- 
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atomic molecules of experimental interests in Refs. I47H57I So far the only nuclear 
AM which has been measured is the AM of the 133 Cs nucleus. It is done by com- 
paring PNC amplitudes between different hyperfine structure sublevels in the same 
PNC experiment where the Cs weak charge is measurecP (the method was proposed 
irPS). Interpretation of the measurements^ indicates some problems. For example, 
the value of Cs AM is inconsistent with the limit on the AM of TlP^ 

To resolve the inconsistencies and obtain valuable information about P-odd 
nuclear forces it would be very important to measure anapole moments for other 
nuclei. In particular, it is important to measure AM for a nucleus with an unpaired 
neutron (Cs and Tl have unpaired protons) . AM of such nucleus depend on different 
combination of the weak interaction constants providing important cross-check. 
Good candidates for such measurements include odd isotopes of Ra, Dy, Pb, Ba, La, 
Lu and Yb. The Ra atom has an extra advantage because of a strong enhancement 
of the spin-dependent PNC effect in the 1 So - 3 ^>2 transition due to proximity of 
the opposite-parity state 3 P? (AE = 5 cm -1 )!^ 

Experimental work is in progress for Rb and Fr at TRIUMF^K] anc j f or rjy 
and Yb at BerkeleyP^E^ Measurements for Xe and Hg are planned at university of 
CreteP 

Calculations of the nuclear spin-dependent PNC amplitudes (including the effect 
of AM) for a number of atoms and ions of experimental interest were reported in 
RefsEDHTDl 

5. Electric dipole moment 

Permanent electric dipole moment (EDM) of neutron, atom or molecule would vio- 
late both P and T invariance. Under conditions of the CPT-theorem this would also 
mean CP- violation. The Kobayashi-Maskawa mechanism of the standard model 
leads to extremely small values of the EDMs of the particles. It is also too weak to 
explain the matter-antimatter asymmetry of the Universe. On the other hand, most 
of the popular extensions to the standard model predict much larger EDMs which 
are within experimental reach. The EDM of an atom or a molecule is mostly due 
to either electron EDM and T, P-odd electron-nucleon interactions in paramagnetic 
systems (with non-zero total electron angular momentum J) or to the T, P-odd 
nuclear forces in diamagnetic systems ( J = 0; nuclear-spin-dependent e-N interac- 
tion contributes here too) . The existence of the T, P-odd nuclear forces leads to the 
T, P-odd nuclear moments in the expansion of the nuclear potential in powers of 
the distance R from the center of the nucleus. The lowest-order term in the expan- 
sion, the nuclear EDM, is unobservable in neutral atoms due to the total screening 
of the external electric field by atomic electrons. It might be possible however to 
observe the nuclear EDM in ions (see below). The first non- vanishing electrostatic 
terms which survives the screening in neutral systems are the octupole moment and 
the Schiff moment. The octupole moment does not produce EDM in diamagnetic 
atoms. Below we discuss the effects of the nuclear and electron EDM and the Schiff 
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5.1. Nuclear EDM 

It was widely believed that one needs neutral particles (e.g., neutron, neutral atom 
or molecule) to study EDMs. This is because the EDM is expected to be very small 
and it would be very hard to see the effect of its interaction with external electric 
field on the background of the much stronger interaction with the electric charge. 
On the other hand, the EDM of neutral systems is very much suppressed by the 
effect of screening of an external electric field by electrons (Schiff theorem - see 
section 5.3). The Schiff theorem may be violated by the relativistic effect (which 
dominates in the case of the electron EDM), the hyperfine interaction and the 
nuclear finite size effect.^ For example, the lowest-order T, P-odd nuclear moment, 
the nuclear EDM is practically unobservable in the neutral systems (except for a 
small contribution due to the hyperfine interaction calculated in Ref. IT2"|) . First 
observable T, P -odd nuclear moment, the Schiff moment, is non-zero due to the 
finite nuclear size. 

It is important therefore to explore the possibility of studying EDMs of charged 
particles (e.g. muons or atomic ions). There are realistic suggestions of this kind 
in Refs. 73(76] based on the use of ion storage rings. There are also proposals to 
measure EDM in molecular ions.^ 

The external electric field is not totally screened on the nucleus of an ion. Its 
value is 

E N = ^E Q , (13) 

where Eo is the external electric field, En is the electric field at the nucleus, Ze 
is the nuclear charge, i^e is the charge of the ion, e is proton charge. The formula 
(|13p can be obtained in a very simple way. The second Newton law for the ion and 
its nucleus in the electric filed reads 

Midi — ZieEo, Mj^aN — ZeEpj, 

where Mi is the ion's mass, is its acceleration, Mjv is nuclear mass (Mjv ~ 
Mj), and cin is its acceleration. Since the ion and its nucleus move together, the 
accelerations must be the same (at — oat), therefore 

E N = ^E ^L « ^E . (14) 
N Z Mi Z V ' 

Quantum mechanical derivations of this formula can be found in Refs. [751501 Nu- 
merical calculations of the screened electric field inside an atomic ion were per- 
formed in a number of our works (see, e.g. Ref. [75| . 

The Hamiltonian of the nuclear EDM (djy) interaction with the electric field is 
given by 

Z- 

Hd = (1nE n = cLn-^Eq. (15) 
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Screening is stronger for diatomic molecules where we have an additional suppres- 
sion factor in eq. <jT4j) , M N /Mi = Mi /(Mi + M 2 ), where Mi and M 2 are the 
masses of the first and second nucleus. For the average electric field acting on the 
first nucleus we obtain 

Zi Mi , . 

^ = ^mTm^°- (16) 

A Quantum mechanical derivation of this formula can be found in Ref. 1801 Note 
that the screening factor here contains both nuclear masses. This indicates that the 
nuclear motion can not be ignored and the screening problem is more complicated 
than in atoms. For example, in a naive ionic model of a neutral polar molecule 
A + B~ , both ions A + and B~ should be located in the area of zero (totally screened) 
electric field since their average acceleration is zero. This could make A + and B~ 
EDM unobservable even if they are produced by the nuclear Schiff moment or 
electron EDM. In a more realistic molecular calculations the Schiff moment and 
electron EDM effects are not zero, however, they may be significantly suppressed 
(in comparison with a naive estimate of ionic EDM in a very strong field of another 
ion) and the results of the calculations may be unstable. 

Note also that contrary to atomic ions the contribution of the Schiff moment to 
the T,P-odd effects in heavy molecular ions exceeds the contribution of the nuclear 

edmP 

In the case of monochromatic external electric field its frequency can be chosen 
to be in resonance with the atomic electron excitation energy. Then for the effective 
Hamiltonian we can have 

H d = d N E N (t) > d N E . (17) 



5.2. Electron EDM 

Paramagnetic atoms and molecules which have an unpaired electron are most sen- 
sitive to the electron EDM. The EDM of such systems can be expressed in the 
form 

d = Kd e , (18) 

where d is the EDM of an atom or molecule, d e is electron EDM, and K is electron 
structure factor which comes from atomic calculations. The factor K increases with 
nuclear charge Z as Z 3 ^ times large relativistic factor R(Za)^ which may exceed 
the value of 3 in heavy atoms. A rough estimate of the enhancement factor in heavy 
atoms with external s 1/2 or p 1/2 electron is K ~ 3Z 3 a 2 R(Za) ~ 10 2 - l0 3 .EHIS2l 

Several orders of magnitude larger K ~ 10 7 — 10 11 exist in molecules due to 
the mixing of the close rotational levels of opposite parity (including A-doublets) 
Following Sandars this enhancement factor is usually presented as a ratio of a 
very large internal molecular field to the external electric field which polarizes the 
molecule. 
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The best current limit on electron EDM comes from the measurements of the 
thallium ED]VpS and reads 

d e = (6.9 ± 7.4) x 10~ 28 e cm. (19) 

Here the value K = — 58!!^ was used for the interpretation of the measurements. 
The value of K for Tl is very sensitive to the inter-electron correlations but three 
most complete calculations^^!] gj ve very close results (there are also calculations 
of Ref. [87] which give smaller result). 

In contrast to paramagnetic atoms the diamagnetic (closed shell) atoms are 
much less sensitive to electron EDM. This is because the only possible direction 
of the atomic EDM in this case is along nuclear spin and hyperfmc structure in- 
teraction must be involved to link electron EDM to nuclear spin. For example, for 
the mercury atom K ~ 10 2 tHHlIHs] However, due to very strong constraints on the 
mercury EDM^l the limit on electron EDM extracted from these measurements is 
competitive to the Tl result (fT9]p Ql 

|d e | < 3 x 10~ 27 e cm. (20) 

New strong limits on the electron EDM were recently found from the measure- 
ments for the YbF molecule!^ Experiments are in progress to measure the electron 
EDM in Cs^Fr^ThO^PbC^and in solid-state experimentsP 

Molecular EDM due to electron EDM were calculated for a number of diatomic 
molecules in Refs.97l fT05l 



5.3. Schiff moment 

Schiff moment is the lowest-order T, P-odd nuclear moment which appears in the 
expansion of the nuclear potential when screening of the electric field by atomic 
electrons is taken into account. This potential can be written as (sec different deriva- 
tions, e.g. in Ref. HM MH HB M 

p(r) 



R 



|R-r| Z 



rd 6 r, 



(21) 



where p(r) is the nuclear charge density normalized to Z, and d is the nuclear EDM. 
The second term in (121[) is the screening. Taking into account the finite nuclear size 
the dipole term in the multipole expansion of (|21l) can be written as^SI 



</>(R) = - 



3S • R 



B 



where B = J p(R)R 4 dr and 



10 



(r 2 v) 



(22) 



(23) 



is the Schiff moment. The expression ([22]) has no singularities and can be used in 
relativistic calculations. More accurate expressions, which include corrections for 
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Table 3. EDMs of diamagnetic atoms of experimental interest. 



z 


Atom 


[S/( e fm 3 )] 
XlO 17 e cm 


T) e cm 
xlO -25 


Experiment 


2 


3 He 


8 x 10" 5 


5 x 10" 4 




54 


129 Xe 


0.38 


0.7 


SeattlefSIlAnn Arbo^ 










Princeton^] 


70 


171 Yb 


-1.9 


3 


BangalorepH KyotrJUn] 


80 


199 Hg 


-2.8 


4 


Seattl^ 


86 


223 Rn 


3.3 


3300 


triumfEHI 


88 


225 Ra 


-8.2 


2500 


ArgonneESlKVff^l 


88 


223 Ra 


-8.2 


3400 





the finite nuclear size ~ Z 2 a 2 , were obtained in Refs. 108, 80] 11091 The authors 
also considered a partial screening of the octupole moment The Schiff moment 
is caused by the T, P-odd nuclear forces. The dominant mechanism is believed to be 
the T, P-odd nucleon-nucleon interaction. Another important contribution comes 
from the EDMs of protons and neutrons. 

Schiff moment is the dominant nuclear contribution to the EDM of diamagnetic 
atoms and molecules. The best limit on the EDM of diamagnetic atoms comes from 
the measurements of the EDM of mercury performed in Seattle^ 

K 199 Hg)| < 3.1 x 10- 29 |e|cm. (24) 

Interpretation of the measurements requires atomic and nuclear calculations. 
Atomic calculations link the EDM of the atom to its nuclear Schiff moment. Nuclear 
calculations relate Schiff moment to the parameters of the T, P-odd nuclear interac- 
tions. Summary of atomicQ3Z>ki!2l anc j nuc i eaj n2H]iiMjrrn>] calculations for diamagnetic 
atoms of experimental interest is presented in Table |3l To compare the EDM of dif- 
ferent atoms we present only the results of our nuclear calculations which all were 
performed by the same method. For Hg and Ra there are several recent nuclear 
many-body calculations available (see references in the most recent calculation^^) . 

The dimensionless constant rj characterizes the strength of the P, T-odd nucleon- 
nucleon interaction (in units of the weak interaction Fermi constant G) which is 
to be determined from the EDM measurements. Using (|24"]) and the data from the 
Table one can get 

5( 199 Hg) = (-1.8 ± 4.6 ± 2.7) x 10~ 13 e cm (25) 
and for the T, P-odd neutron-proton interaction 

?7„ p = (1±3±2) x 10 -5 . (26) 

5.4. Nuclear enhancement. 

In some nuclei with the quadrupole deformation there is a close opposite parity 
level with the same angular momentum. This leads to an enhancement of nuclear 
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and Schiff momentP-^ U p to an order of magnitude. Significantly larger 
enhancement is possible in nuclei with the octupole deformation.^^ This can be 
explained in a very simple way. Nuclear deformation creates a collective intrinsic 
Schiff moment in the rotating nuclear reference frame 

Sintr » eZR%-J^L=, (27) 

207Tvo5 



where Rn is the nuclear radius, 02 ~ 0.2 is the parameter of the quadrupole 
deformation, and rs 0.1 is the parameter of the octupole deformation. The 
intrinsic Schiff moment (|27|) does not violate T or P invariance and, if no T, P- 
odd interaction is present, it averages to zero in the laboratory reference frame 
due to the nuclear rotation. However, when T, P-odd interaction is included, it can 
mix close rotational states of opposite parity (which are similar to f2- doublet in 
diatomic molecules). Small energy interval between these states leads to a strong 
enhancement of the nuclear Schiff moment in the laboratory reference frame 

Siat ~ ( t |gPT i~ ) S intr ~ 0MeP 2 piZA 2 ISrl =-^= ~ 700 x 10~V fm 3 , 

_c/-}_ — Hi— hi+ — hi— 

where ro = 1.2 fm is the inter- nucleon distance \E + — E-\ ~ 50 keV. The estimate 
([28)) is about 500 times larger than the Schiff moment of a spherical nucleus like 
Hg. 

It was pointed out in Ref. II 181 that the octupole deformation doesn't need to be 
static. Soft octupole vibrations lead to similar enhancement. 

Large values of the Schiff moment for Ra and Rn (see Table [3]) are due to the 
nuclear octupole deformation. Some additional enhancement of the atomic EDM 
is due to the large nuclear charge. The nuclear enhancement may also manifest 
itself in molecules, e.g. in RaOp^ where the T, P-odd nuclear spin- molecular axis 
interaction exceeds that in the experimentally studied TlpESElH about 500 times. 
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